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CD ' We review lattice evidence showing that center-vortex condensation is a serious 

Q\ ' candidate for the mechanism of colour confinement in quantum chromodynamics. 

o' 

o 

f^ ' 1 Introduction 

j^ , Any model that contends for the true theory of colour confinement in SU(-/V) 

gauge theory has to be able to provide an explanation of or an answer to the 
^ I following facts and questions: 

• All non-zero iV-ality colour charges are confined. 

• All zero iV-ality charges are screened at large distances. 

• At intermediate distances, potentials between higher representation char- 
^1^ ■ ges rise approximately linearly, and the corresponding string tensions are 

H \ proportional to the eigenvalues of the quadratic Casimir operator of the 

representation. (We dubbed this phenomenon Casimir scaling.) 

• A deconfinement transition occurs at high temperatures or high densities 
of hadronic matter. What is the nature of the transition? What hap- 
pens to confining configurations at the transition, what symmetry gets 
restored or broken? 

• What is the relation to other phenomena, like chiral symmetry breaking? 
In this talk I will argue that confining configurations that have a chance to 
successfully address the above issues, are center vortices. I will mostly review 
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recent results of our group obtained in lattice investigations of the confinement 
problem. Other aspects of the vortex picture were covered at this Conference 
in the plenary talk of ReinhardtEl (see also short contributions of AlexandruH, 
Engelhardta, and Langfelda). 

The vortex-condensation model c£ confinement stems from late seventies, 
and was formulated first by 't HooftQ and developed by many authors. Let 
me briefly recall the original idea. To characterize phases of the gauge theory, 
Wilson Q suggested to use the loop operator 



A{C) = Tr 



P exp I ig (p Apdx^ j 



(1) 



In a pure gauge theory, confinement is signalled by the leading behaviour of 
this quantity, namely 

TTr/^\ I A/^\\ { e^'^^^'~^> ... confinement, .„s 

^^^) = ^^^^))-[e-.nc, ... deconfinement, ^^^ 

where .4(C) is the area of the minimal surface spanned by the loop, V{C) its 
perimeter. Wilson's loop operator measures the magnetic fiux through C and 
creates electvic fiux along C. 

't HooftEl suggested another operator, measuring the electric fiux through 
some other loop C" and creating magnetic fiux along C". It was defined through 
commutation relation with A{C) {n is the linking number between C and C): 

A(C)B{C') ^ B(C')A{C) exp (27rm/iV) . (3) 

The expected behaviour of 't Hooft's operator is opposite to the Wilson loop: 

/R(r^'\\ /e"''^(<^') ... confinement, , . 

{^(^))-<^^^,Aic") ... deconfinement. ^^ 

The effect of the 't Hooft loop operator B(C') on physical fields is a gauge 
transformation which is singular on the curve C . It creates a thin vortex of 
magnetic flux (line-like in 3D, surface-like in 4D). To avoid infinite energy, the 
singularity along the loop C" has to be smoothened: vortices acquire a core of 
certain thickness. Many things in this talk will be related in some way to the 
thickness of the vortex core. 

The essence of the vortex model of confinement is that the QCD vacuum is 
filled with closed magnetic vortices that have the topology of tubes (in 3D) or 
surfaces (in 4D). Thick vortices condense in the vacuum. A simple argument 
can explain the area law for the Wilson loop through fiuctuations in the number 
of vortices linked to the loop. 



It is a sort of paradox that this idea has for years not been subjected to 
extensive numerical tests in lattice simulations. However, in the last couple of 
years the idea returned to the stage and got into the focus of attention, 
turn was due to the discovery of center dominance in maximal centeiigaiigeC 
and the following accumulation of evidence in favour of the pictureEl lj. 

The outline of my talk is the following: In Section I will briefly sketch 
our recent strong-coupling calculationEj, which shows that center vortices are 
stable saddlepoint configurations of a long-range effective action that can be 
derived from the usual lattice Wilson action. Then (Section ||), I will introduce 
our method to detect center vortices that is based on center projection in max- 
imal center gauge, and review some interesting results obtained in numerical 
simulations. Section |4| then addresses a question, raised recently by Bornyakov 
et al.Ej, how physical results depend on the number of gauge copies used in it- 
erative MCG fixing. Finally, I will argue in Section g that the observed Casimir 
scaling can be accommodated with the vortex picture, and that the thickness 
of the vortex core is vital for the explanation of the scaling behaviour. 

2 Confining Configurations of the Long- Range EfTective Action 

There are various ways to argue for the center of the gauge group being impor- 
tant for the confinement mechanism. The asymptotic string tension for charges 
in the representation r depends only on its iV-ality (i.e. Z^ charge); the decon- 
finement transition is associated with spontaneous breaking of the global Z^ 
symmetry; etc. We have recently provided another argument: if one calculates 
an effective action at larger distances (and strong coupling) starting from the 
usual Wilson action, its stable saddlepoints are thick center vortices. 

The calculation proceeds in the following wayllj: We start from the pure 
gauge theory at strong coupling on a fine lattice (links on the fine lattice are 
denoted generically as U) , and try to calculate the effective action on a coarse 
lattice with links V ^ where the coarse lattice spacing is L times the original 
L'^-lattice spacing. The idea is to derive an effective action where the leading 
contributions to any Wilson loop on the 1^-lattice are obtained from a local 
action. To achieve this, we integrate over all links on the [/-lattice except 
[/-links on small 2-cubes around sites of the V^-lattice, as shown in Fig. Il|.'^ 

This defines effective action Sl'- 



Z= f[DV] f n idUi]\ f n [dUr]l[s[v^]{UU..U),~I 

2-cuboEi ^_ 2-cubcs 

"We work in 3 dimensions. The extension to 4D should be straightforward. 
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Figure 1: Unintegrated links on the coarse V-lattice. 



/[^^]/ n [dUi]exp(^SL[V,U] 



/e2-cubes 



(5) 



Instead of. Ui we introduce in the 2-cubes group- valued plaquette variables 
/i, g (Fig. 0) EZl. After integrating out the g variables, the result is (displaying 
only terms of low order in both h and /?): 



Z « J[DV] [Dh] n { 1 + 2 (0 ' E Xi [(hhh) 



-'! 



X exp 



fETrW + 2ff 



E Xi[(hhh)cAhhh)c 
adj. ciC2eK 

I' 

r 2 

+2(1) T.^r[VVV^V^] 
\ / p, 



(6) 



This resembles an adjoint-Higgs Lagrangian, with an SU(2) gauge field 
V^ coupled to 24 "matter" fields h in the adjoint representation. Note that 



2-Cube 



h,, 





Figure 2: Plaquette variables on the 2-cubes. 
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for large L, the "Higgs" potential term is much larger than the "kinetic" and 
pure-gauge (y-plaquette) terms, so the /i-fields fluctuate almost independent 

We then make a unitary gauge fixing of the /i-fields (breaking SU(2) to 
Z2), integrate out the remaining h d.o.f., and arrive at: 

4(L-2) 



2(f) E//'''Tr[(4).T^'(/.L)'.^I 



+ 2(7) Et^[^^^^^^]- (7) 

\ / p, 
Let us now look for saddlepoints. Sunk is maximized by 

V^{n) = Z^{n) X 5(n).gt(n + Ai), Z^ = ±1, (8) 

where g{n)g'' {n + /j,) is fixed by the particular unitary gauge choice, while Spi^q 
is maximized if ZZZZ = +1. This is the unitary gauge ground state. 

One can create a thin center vortex on this state by a discontinuous gauge 
transformation, e.g. 

Zyin)-[ll Ttherir^'' Mn)^Z.in) = l. (9) 

This configuration is stationary: Sn^k is still a maximum (insensitive to cen- 
ter), Spiaq is extremal (maximal or minimal) on all plaquettes. 
Stability depends on the eigenvalues of 



^ ^eff _ S"^ Slink , ^ ^plaq 



(10) 



SV^{ni)SVAn2) SV^ini)dV,{n2) SV^{ni)SVAn2)' 
where 

s'sunk /^y(^-^)+^^ 6^Spi,q f^Y^ 

SVSV V4/ ' SV6V U/ ■ 

Crucial fact is now the following: for /3/4 <C 1 and 

4:{L - 2) + 12 < L^ =^ L>5 (12) 

the contribution of S'^Spiaq/SVSV to the stability matrix (and therefore to its 
eigenvalues) is negligible compared to S'^Sunk/^VSV. 
This has (at least) three important implications: 



• Vortex stability: The thin vortex is a stable saddlepoint of the full effec- 
tive action S^jj at L > 5. 

• Vortex thickness: A "thin" vortex on the ^-lattice means thickness < L 
on the [/-lattice. This means that stable center vortices are « 4 — 5 lattice 
spacings thick. This is the distance where the adjoint string breaks at 
strong couplings! 

• Percolation: From S^ff, we see that center vortices in _D = 3 cost action 

8 {(3/4:) /unit length, while the entropy is 0(l)/unit length. 

Entropy ^ action implies that vortices percolate through the lattice, and 

confine nonzero A^-ality charges. 
However, the above result has been obtained in the strong-coupling ap- 
proximation. To study the role of vortices at weak couplings, we resort to 
numerical simulations. 

3 Center Projection in Maximal Center Gauge 

The model of 't Hooft relies on no particular gauge. However, without gauge 
fixing, one can derive certain relations that, at the fip.t sight, seem related 
to vortices, but apparently have no physical relevancell3. Gauge fixing turns 
out to be also useful in pin-pointing the relevant (for vortices and confine- 
ment) degrees of freedom. One can then study in detail e.g. their topological 
propertiesll3, behaviour across deconfinement phase transitiontHi, etc. 

We proposed tflj-identify vortices in thcrmalized lattice configurations in 
the following steps QB: 

First we fix to maximal center gauged by maximizing the expression: 



• 



7^^^|Tr[^7^(a;)]| . (13) 

This in fact is adjoint Landau gauge; the above condition is equivalent 
to maximizing 

7e'^^Tr[t/^^(x)]. (14) 

• Then we make center projection by replacing: 

U^,{x) ^ Z^{x) = sign Tr[t/^(a:)] . (15) 



• 



Finally we identify excitations (P-vortices) of the resulting Z2 lattice 
configurations. 



''For the continuum counterpart of MCG see the talk of Reinhardt'^'^^ 



A series of findings shows that vortices identified via center projection in 
MCG, are crucial for colour confinement and related phenomena: 

• Center dominance: It was observed both, in SU(2) ErH and (less con- 
vincingly) in SU(3) lattice gauge theory E3 that the string tension ob- 
tainecLfrom center-projected configurations in MCG (or Laplacian center 
gauge L^) agrees remarkably well with the asymptotic string tension of 
the full theory. 

• P-vortices locate center vortices: u One can define "vortex- limited" Wil- 
son loops Wn{C), i.e. Wilson loops evaluated on the original, unpro- 
jected lattice, on a subensemble of configurations in which n P-vortices 
pierce the minimal area of the loop. One expects Wn{C)/Wo{C) to 
approach (—1)" asymptotically; this expectation is confirmed for large 
enough loops in numerical simulations. 

• P-vortices locate physical o&jects;Ej'u Vortex density scales according to 
the asymptotic freedom, like a dimensionful quantity ~ A^. 

• Center vortices are correlated with confinement: D'lI Removal of center 
vortices destroys confinement. 

• Center vortices are correlated with chiral symmetry breaking: u Removal 
of center vortices restores chiral symmetryj.^ 

• Center vortices cause non-trivial topology :u Lattice configurations with- 
out vortices belong to the trivial topological sector (have zero topological 
charge) . 

• Decoufiwment can he understood as a center vortex percolation transi- 
tion: Il3't3 In the confined phase, most P-plaquettes belong to a single 
huge vortex cluster percolating through the whole lattice volume. In the 
deconfined phase, vortex clusters cease percolating in space slices, while 
they continue percolating in time slices. Most vortices are short loops, 
winding in the time direction. In this way, perimeter-law behaviour of 
time-like Wilson loops can peacefully coexist with area-law behaviour of 
space-like Wilson loops. 

4 Center Dominance, Gauge Copies, and Lattice Size 

Though the above list of indications in favour of our vortex identification pro- 
cedure looks quite impressive, a word of caution is necessary. MCG, similar 
to Landau or maximal abelian gauge, suffers from the Gribov copy problem. 
The iterative procedure, used for gauge fixing, converges to a local maximum 
which is slightly different for every gauge copy of a given lattice configuration. 
At the first sight, Gribov copies in MCG do not seem to be a severe problem 
in our procedure; it appears that P-vortex locations vary comparatively little, 
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Figure 3: The dependence of the Creutz ratio x(4, 4) on the number of gauge copies Ncopy 
for /3 = 2.5, for various lattice sizes. 



from copy to copyB. 

However, a serious drawback might be a strong dependence of physical 
results (Creutz ratios, vortex density, etc.) on the number of gauge copies 
used in MCG maximization procedure. For each lattice configuration one can 
generate a set of its random copies, fix all these copies to MCG, and evaluate 
quantities of interest on the "best" copy, that is the one with the highest 
value of TZ in Eq_,(|l3[). Such an investigation has recently been reported by 
Bornyakov et al.Ej They e.g. show that the Ncopy —* oo extrapolation of the 
center-projected string tension underestimates the full string tension by as 
much as 30% (at /3 = 2.5)! This seems to indicate that the observed center 
dominance in MCG (obtained from a small number of copies) might have been 
just a numerical coincidence. 

In our opinion, the reported strongjiisagreement is due to finite size effects. 
The lattices used by Bornyakov et al.Ej seem rather small compared to those 
used in our earlier studies. In fact, what is "rather small" or "large enough" 
depends on what are the typicaL|Sizcs of vortex cores. Independent sources 
of estimating vortex thickness la'c2lE3 point towards a value of a little over one 
fermi, i.e. center vortices are w 12-14 lattice spacings thick at /3 = 2.5 and 
the lattice size, used by Bornyakov et al. (16^), seems too small. We have 
therefore repeated their analysis on lattices of various sizes; it turns out that 
on sufficiently large lattices the dependence on the number of gauge copies is 
much weaker, and center dominance is quite accurate. This fact-is illustrated 
in Figs. and H, details can be found in our recent publicationcJ.R 

'^This might not be the final answer to the problem. We have recently been informed that, 
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Figure 4: Center-projected Creutz ratios Xcp{I, I) vs. /, from our largest lattices. 



5 Casimir Scaling 

We have repeatedly emphasized that trying to understand the confineinent 
mechanism one has to be able to explain the problem of Casimir scaling. If 
we look at string tensions for higher representation colour charges, we find 
different behaviour in two distinguishable regimes: 

• Casimir- scaling regime: At intermediate distances, from the onset of 
confinement to the onset of screening, the string tension is roughly pro- 
portional the quadratic Casimir. For SU(2), e.g., 



0'l/2 



-,3{J + !)■ 



(16) 



• Asymptotic regime: The SU(iV) quark colour charge is eventually screen- 
ed by gluons to the lowest representation with the same transformation 



using a combination of simulated annealing and the usual (over)relaxation, one can still find 
copies with higher maxima of TZ [V. Bornyakov, M. Polikarpov, private communication]. 
Also, using[tp.ndau gauge preconditioning, one reaches a higher maximum, but loses center 
dominance Ej- Though we were [*Hple to find a cause of this problem, the loss of what we 
called "vortex-finding property" c3 during LG fixing, it is clear that a better gauge-fixing 
procedure and/or a imndjftcation of the gauge-fixing condition, is needed. Some alternatives 
have been proposed EJ'ESq. 
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properties under Zjv- For SU(2): 

(7i/2 • ■ • j = half-integer, 



... J = integer. ^ ' 

As N —>■ cxD, colour screening is suppressed, Casimir scaling is_exact. But even 
at A^ = 2, 3 there is a Casimir-scaling region of finite extent E3Ej. 

Center vortices are good at explaining the asymptotic region. Loops which 
transform trivially under the group center (such as the adjoint) are unaffected 
by fluctuations in the center elements, and therefore shouldn't get an area law. 
This is strictly true, however, only for thin center vortices (P-vortices) . 

Is there a way of reconciling approximate Casimir scaling at intermediate 
distances with center vortices? The answer is apparently yes, and the crucial 
ingredient in the explanation of Casimir scaling through vortices is the vortex 
thickness. If the slice of a thick vortex, in the plane of a Wilson loop, is entirely 
inside (or entirely outside) the loop, the effect is the same as in the thin vortex 
model. However, if the vortex overlaps the loop perimeter, one can envisage its 
effect on the loop by insertion of a group element G, interpolating smoothly 
between —1 and 1 (as would be the case in an abelian theory). Were the 
vortices relatively thin on average, this effect would be unimportant; it would 
lead to perimeter-like correction to Wilson loops up to sizes comparable to 
the size of the vortex core. But, as already discussed in the previous section, 
there are good reasons to believe that vortices are quite thick. This enables 
one to explain the gross features of the behaviour of string tensions both in 
the Casimir-scaling ani-asymptotic-regions. A simple model has been worked 
out in oui:, earlier workEil (see alsoE3). 

Balic3 and Deldarcj have recently published results of high statistics com- 
putations of higher representation potentials in the SU(3) lattice gauge theory. 
Casimir scaling turns out to be quite precise, much more precise than in our 
over-simplified model. We believe this is not a real problem for center vortices, 
and our model and Ansatz for vortex-core profiles can be refined to explain 
the observed behaviour. 

Finally I want to mention that to be compatible with N —> oo where 
Casimir scaling is exact, the vortex thickness should grow with M.. Some 
evidence for such a growth has recently been presented by Monterol23. 

6 Epilogue 

The first in this series of Conferences took place in 1994 at the Lake of Como. 
There is no mention of the vortex model of confinement in its Proceedings. 
The idea, after its happy infant years, was sleeping deeply at that time like the 
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(a) E. Burne-Jones (1833-1898): The 
Sleeping Princess (Briar Rose), detail. 



(b) J. W. Waterhouse (1849-1917): 
The Lady of Shalott, detail. 



Figure 5: "Confinement I" (1994) versus "Confinement IV" (2000). 

Briar Rose (Fig. Ha). I tried to convince you in my talk that now, six years 
later, the Sleeping Beauty is up after a long sleep, and looks very attractive 
(Fig. ||b). Some problems, however, persist. I sincerely hope that, in the years 
to come, she will not fall asleep again, or turn into the Maleficent Queen. 
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